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Abstract. In this paper, the dynamics of spontaneous shape fluctuations of a single, giant quasi-spherical 
vesicle formed of a single lipid species is revisited theoretically. A coherent physical theory for the dy- 
namics is developed based on a number of fundamental principles and considerations and a systematic 
formulation of the theory is also established. From the systematic theoretical formulation, an analytical 
description of the dynamics of shape fluctuations of quasi-spherical vesicles is derived. In particular, in 
developing the theory we have made a new interpretation of some of the phenomenological constants in a 
canonical continuum description of fluid lipid-bilayer membranes and shown the consequences of this new 
interpretation in terms of the characteristics of the dynamics of vesicle shape fluctuations. Moreover, we 
have used the systematic formulation of our theory as a framework against which we have discussed the 
previously existing theories and their discrepancies. Finally, we have made a systematic prediction about 
the system-dependent characteristics of the relaxation dynamics of shape fluctuations of quasi-spherical 
vesicles with a view of experimental studies of the phenomenon and also discussed, based on our theory, a 
recently published experimental work on the topic. 

PACS. 68.15.+e Liquid thin films. - 87.I6.Dg Membranes, bilayers, and vesicles. 



1 Introduction 

Dynamics of fluid membranes was already observed a cou- 
ple of centuries ago, as "flickering" seen under optical 
microscope of resting red blood cells. The first physical 
study of the "flickering" phenomenon of red blood cells 
was made, however, only in 1975 by Brochard and Lcnnon 
jj]]. They presented the first theory for the phenomenon, 
where they employed the concept of bending elasticity due 
to Helfrich Q for describing the conformational flexibility 
of fluid membranes and associated the "flickering" with 
the dynamics of equilibrium shape fluctuations of red-cell 
membranes that are thermally driven due to the conforma- 
tional flexibility of the membranes. This association has 
since become one of the essential ideas in theories for dy- 
namics of fluid membranes. 

During the last three decades, unilamellar lipid-bilayer 
vesicles reaching micrometer sizes that can be routinely 
prepared in the laboratory have provided one of the most 
important classes of model systems for biological cell mem- 
branes. They allow for well-defined and quantitative phys- 
ical investigations of the conformational behaviour of fluid 
membranes. The relevant physics governing the static equi- 
librium shapes of giant, unilamellar vesicles has been fully 
understood, as the result of a series of both theoretical 
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and experimental studies ||,||,[|,|]]. The static properties 
of equilibrium shape fluctuations have also been investi- 
gated The understanding of the dynamics of equi- 
librium shape fluctuations of the vesicles has not, however, 
reached the same stage. This is largely the consequence of 
the fact that the dynamics of fluid membranes involves 
many more physical factors than the statics. The static 
properties of a fluid membrane in thermodynamic equilib- 
rium are entirely determined by thermodynamic energet- 
ics of the membrane alone. In contrast, the dynamics of 
the membrane is not only related to the thermodynamic 
energetics, but also depends on the dynamics of the bulk 
aqueous fluid in which the membrane is suspended: mo- 
tion of the membrane induces motion in the bulk fluid, 
which in turn exerts hydrodynamic force on the mem- 
brane. Furthermore, any intrinsic dissipation mechanisms 
of the membrane will in principle play a role in the dy- 
namics. 

Many studies, both experimental and theoretical, have 
actually focused on the dynamics of fluid membranes. The 
earlier work of Schneider et al. ]l(| investigated experi- 
mentally by use of fluorescence microscopy the fluctuation 
spectrum of giant, unilamellar and quasi-spherical vesi- 
cles. In the theoretical interpretation of the experimental 
data, Schneider et al., and later, Milner and Safran |Q, 
treated the energetics of the membrane as that of a single 
and incompressible fluid surface with the Helfrich bending 
rigidity. Moreover, they neglected any intrinsic dissipation 
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of the membrane. Their theory predicted a single relax- 
ation rate for each spherical-harmonic mode of the vesicle 
shape fluctuations. 

The earlier notion that a lipid-bilayer membrane could 
be considered as a single and incompressible fluid surface 
soon had to be modified. The relevance of both the bi- 
layer architecture of the membrane and the compressibil- 
ity of the constituting monolayers was first appreciated 
and demonstrated in the context of equilibrium shapes 
of giant lipid-bilayer vesicles ||||,|ll],|l^,[l3|. Evans and 
Yeung soon also pointed out the relevance of this new 
notion to the conformational dynamics of giant vesicles 
EJ . Based on their analysis of a tether-pulling exper- 
iment on giant vesicles, they proposed that the bilayer 
structure of the membrane in fact allows the two mono- 
layers to move relative to each other laterally and that 
the relative motion is associated with dissipation of en- 
ergy phcnomcnologically described as intermonolayer fric- 
tion. Moreover, they suggested a coupling between bend- 
ing of the membrane and a local field of relative area di- 
lation/compression of the monolayers, whose relaxation 
in turn drives the relative intermonolayer motion. In this 
theory, therefore, the conformational dynamics of a fluid 
membrane would be coupled to dynamics of the relative 
monolayer dilation/compression field. 

Following the idea put forward by Evans and Ycung 
PI , Seifert and Langer formulated a theoretical descrip- 
tion of the paradigmatic case of the dynamics of equilib- 
rium fluctuations of an almost planar bilayer membrane 
embedded in a viscous aqueous medium They also 
discussed the relevance of their result to the experimen- 
tal measurements made by spin-echo technique on mul- 
tilamellar membrane stacks. No conclusive confirmation 
of the theory, however, could be made due to the lack of 
direct experimental data. 

Giant, quasi-spherical vesicles, whose fluctuations can 
be observed directly by video microscopy, provide another 
type of systems, for which a theory for the dynamics can in 
principle be put to a quantitative test against experimen- 
tal data. Yeung and Evans |l6| , and later, Bivas et al. |T7j , 
considered such systems. Each group presented a theoret- 
ical formulation of its own for the dynamics of equilibrium 
fluctuations of quasi-spherical vesicles. Although, qualita- 
tively speaking, the two theories were similar in that both 
used the ideas originally proposed by Evans and Yeung 
]T^| , they differed both in their detailed descriptions of 
the thermodynamic energetics of the vesicles and in their 
formulations of the equations of motion. Consequently, 
the two theories yielded different analytical results that 
characterize the dynamics. It is, however, not straight- 
forward to justify either of the two theories. Yeung and 
Evans employed a description of the thermodynamic free 
energy of the vesicle membrane, that would be consistent 
with the assumption that local density inhomogeneities 
in the monolayers relaxed much faster than bending de- 
formations. But the assumption is inconsistent with the 
basic starting point of their theory, which was to treat the 
dynamics of bending and density inhomogeneities on the 
same footing. Bivas et al., while proposing and using a 



phenomenological model for the thermodynamic free en- 
ergy that was consistent with their considerations of dy- 
namics, adopted an ad hoc approach to the formulation of 
equations of motion that ultimately governed the dynam- 
ics. The approach lacks the rigor that a systematic ap- 
proach based on some fundamental principles would oth- 
erwise have. Its justification is, therefore, less than trans- 
parent. Moreover, a number of approximations were used 
in both theories, for which the justifications were not clear 
either, due to the absence of a systematic formalism. It is 
thus our opinion that a physically consistent and system- 
atic theory for the dynamics of equilibrium fluctuations of 
quasi-spherical vesicles is still needed. 

Continuing experiments on giant, quasi-spherical vesi- 
cles by the use of phase contrast video microscopy fur- 
ther underline the need for a well-defined theory. A recent 
analysis [jl8) of the experimental data obtained by the use 
of the techniques indicates that the dynamics of equilib- 
rium fluctuations of giant vesicles involves more relaxation 
processes than the single process of relaxation of bending 
predicted by the theory of Schneider et al. jlO| and Mil- 
ner and Safran §. The still improving time resolution of 
the video-microscopy technique will also lead to a large 
amount of experimental data of better quality. Having a 
systematic and well-defined theory will greatly facilitate 
the analysis and interpretation of the experimental data 
as well as the quantitative test of the theory itself, thus 
furthering our understanding of the phenomenon. 

The purpose of this paper is, therefore, two-fold. First, 
we will present, based on a number of fundamental prin- 
ciples and considerations, a coherent physical theory, and 
a systematic formulation of the theory for the dynamics 
of equilibrium fluctuations of giant, quasi-spherical vesi- 
cles formed of a single lipid species. One fundamental 
consideration underlying the theory is Onsager's regres- 
sion hypothesis: the average regression of the spontaneous 
thermal fluctuations in a macroscopic system obeys the 
sames laws as the corresponding macroscopic irreversible 
processes Jl9|| . Following this hypothesis, we may under- 
stand the dynamics of the thermal fluctuations in a gi- 
ant vesicle by investigating the near-equilibrium dynam- 
ics of the same system. Our starting point for formulating 
the near-equilibrium dynamics of a giant vesicle is then 
to describe the vesicle dynamics as the hydrodynamics 
of two coupled systems of (quasi-)two-dimensional fluid 
with time-dependent surface geometry. This description 
of the surface hydrodynamics of the membrane thus nat- 
urally involves both a thermodynamic description of the 
membrane and a description of the hydrodynamics of the 
bulk fluids surrounding the vesicle from within and with- 
out. Formulated this way, the theory, including the vari- 
ous approximations that will be made, appears conceptu- 
ally transparent and systematic. This feature not only will 
make the comparison between the theory and the exper- 
iments, and in turn, any necessary modifications easier, 
but will also lend the formalism readily to any extension 
that deals with vesicle systems of more complexity, for 
example, giant vesicles made either of two lipid species or 
of a lipid species with a membrane protein species. Sec- 
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ondly, we will present, as an application of the general 
formalism, a specific derivation of analytical expressions 
for the regression dynamics of shape fluctuations of quasi- 
spherical vesicles. In particular we discuss systematically 
the different approximations employed in the application. 
Our results are different from those derived by Yeung and 
Evans jl6| and by Bivas et al. |l7j , and we will discuss the 
differences. To facilitate the comparison of the theory with 
the experiments, we also present a systematic numerical 
analysis of the analytical results. 

The plan of the presentation is as follows. In Section 
2, we introduce and discuss a thermodynamic description 
of fluid lipid-bilayer membranes. In Section 3, we present 
a systematic formulation of near-equilibrium dynamics of 
fluid-bilaycr vesicles, which can in principle be applied to 
vesicles with arbitrary equilibrium shapes. In Section 4, 
we apply the general formulation described in Section 3 
to systems of quasi-spherical vesicles. Finally, in Section 
5, we discuss the implications that our results have for ex- 
perimental measurements on the dynamics of spontaneous 
shape fluctuations of quasi-spherical vesicles. We will also 
discuss the differences and connections between our re- 
sults and the theoretical results obtained by other groups 
previously p||l7|l. 

2 Thermodynamic description of fluid 
lipid-bilayer membranes 

In this section, we will consider a vesicle formed of a closed 
fluid lipid-bilayer membrane containing only a single lipid 
species. The vesicle encloses a volume of V, which can 
be considered fixed under typical experimental conditions. 
The two constituting monolayers of the vesicle membrane, 
the outer one labeled by superscript "+" and the inner 
one by "— " , contain N + and N~ lipid molecules, respec- 
tively. As has been recognized ||, N + and N~ depend on 
the vesicle formation process and may be different from 
each other. A total can be resolved with 

techniques of microscopic resolutions or with equivalent 
means, can be defined for the whole membrane. To a good 
approximation, A m may be treated, for most of the mem- 
brane systems studied, as a fixed quantity that depends 
on N + and N~ . When a vesicle is flaccid, in other words, 
having a non-zero excess area A = A m /4irR 2 — 1, where 
R v = (3V r /47r) 1 / 3 , it undergoes thermally driven confor- 
mational fluctuations as a result of the membrane flexibil- 
ity. In the case of giant vesicles that are typically investi- 
gated, whose sizes are in the micrometer range and whose 
bending rigidities are of the order of 10 -19 J, the fluctua- 
tions remain controlled. Thus, the vesicle we consider has 
a well-defined equilibrium shape, whose specific geometry 
depends on V, N + , N~ and temperature T. 

We will present a phenomenological description of the 
thermodynamics of a fluid lipid-bilayer membrane, which 
is consistent with our starting point that the membrane is 
conformationally flexible and consists of two weakly cou- 
pled compressible monolayers. This description will form 
a necessary part of the formulation of the membrane dy- 
namics. 



2.1 Notations and relevant fields 

To describe the conformational flexibility of a fluid mem- 
brane, we use the basic notation of surface differential ge- 
ometry. At length scales larger than molecular lengths, 
the varying conformation of a fluid membrane suspended 
in a bulk aqueous fluid can be approximated by a two- 
dimensional surface embedded in the three-dimensional 
space and with a time-dependent geometry. Such a sur- 
face can be represented by a three-dimensional vector, 
-R^ 1 , it 2 ,t), where u 1 and u 2 are the two coordinates 
parametrizing the two-dimensional internal space and t 
represents time. The derivatives of R(u l , it 2 ,t) with re- 
spect to the internal coordinates u a (a = 1,2) define two 
local tangent vectors to the surface: 

ta (u\u\ t) ee dR (<f^ = da R{u\u\ t), a =1,2. 
ou a 

w 

From t\ and £2, a local unit vector normal to the surface 
can be constructed: 

n = ; r ■ 2 

|*i x t 2 | 1 ' 

The metric tensor a a/ 3 of the surface is then given by 

a af3 = tot * tp J (3) 

which is a symmetric tensor with a non-zero determinant 
a = det(a Q| a). A local area element, dA, can be expressed 
as 

dA= ^fadi^du 2 . (4) 
The inverse of the metric tensor, a*'', follows from 

a Q % 7 = 8° , (5) 

leading to the definitions of contravariant surface tangent 
vectors t a = a a ^tp. The rule of summation over repeated 
Greek indices has been used in the above expressions, and 
will be kept to in the rest of the paper. Clearly, t a ■ tp — 
5p. Any space vector w tangent to the surface can be 
expressed as w = w a t a = w a t a , where w a — w ■ t a and 
w a = w ■ t a = a a ^wp are the covariant and contravariant 
components of w, respectively. 

A full description of the essential local properties of 
the surface requires also the information contained in 

d a d R = r^pd^R + b a(} n , (6) 

where b a p = n-d a dpR is the so-called "curvature tensor", 
and defines the Christoffel symbol. From the curva- 
ture tensor both the mean curvature H and the Gaussian 
curvature K can be obtained, respectively, 

and 

K ee det (a"^) = i- ■ i- (8) 
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where R\ and R2 are the two principal radii of curvature. 

By use of the Christoffel symbol J^L , covariant differ- 
entiations with respect to u a , denoted by D ai of surface 
scalars, vectors and tensors can be performed. For any 
surface scalar, p(u 1 ,u 2 ), Dap^ 1 ,u 2 ) — daplu 1 ,u 2 ), and 
for a surface vector represented by its contravariant com- 
ponents, W a , DaW 1 = daW 1 + r^pW 13 . 

In order to fully characterize the physical state of a 
flexible and compressible fluid lipid-bilayer membrane formed 
of a single lipid species, two more fields are needed in ad- 
dition to i^it 1 , it 2 ,t): p + (u 1 ,u 2 ,t) and p~ (it 1 , u 2 , t) , rep- 
resenting the local surface (number) densities of the outer 
and inner monolayer, respectively, and defined with re- 
spect to the surface described by R^ 1 ,u 2 ,t). 



into account a possible coupling between the bending of 
the membrane and the difference between local monolayer 
surface densities: a non-zero (p + — p~) may be viewed as 



amounting to a dynamic spontaneous curvature 
Eq.fl9h can be rewritten as 



L5L2C 



F 



dA 



2.2 Thermodynamic free energy 

Imagine a vesicle suspended in a bulk fluid at tempera- 
ture T. Its equilibrium (relaxed) state is then character- 
ized both by an average shape Rq{v},u 2 ) and by average 
monolayer densities pj (u 1 ,^ 2 ) = N + /A and (u 1 , u 2 ) = 
N~ /Aq [], where A is the area denned by the average 
shape. For later use, we define 



N+ + N- 



N A = 



N+ - TV- 



Assume now that the vesicle is in a near-equilibrium state 
with an effective shape described by i?(u 1 , u 2 , t) and with 
density fields p + (u 1 ,u 2 ,t) and p~ (u 1 ,u 2 ,t) describing re- 
spective inhomogeneities in the two monolayers. We as- 
sociate the following thermodynamic free energy to the 
near-equilibrium state of the vesicle: 



F 



dA 



P --l 

Po 



dA ^L f ^ 1 



P_ 

Po 



+ / dA 7 (p+ -p )+ / dA 7 (p - po) + 



dA 



dA ^(2H) 2 



dA Atf(— - — ) 



Po 



Po 



(9) 



The five terms in the first two lines may be viewed as 
the result of an expansion of a local free energy function, 
Ft,a(N%,Na,Aa) — AaJ(T, p+ ,p~), about a particular 
reference value, po, of the monolayer surface densities. In 
other words, the first four terms describe the free-energy 
cost associated with deviations in the monolayer surface 
densities from the reference value and the last of the five 
terms represents the free energy required to change the 
area of the effective shape of the vesicle at constant den- 
sity, po- Consequently, the constant parameters, k c s, 7, 
and a, all depend on the chosen reference state. k e g, in 
particular, is the monolayer compressibility modulus cor- 
responding to the chosen po. The first term in the last 
line of Eq.M represents the free-energy cost associated 
with bending of the membrane, and the second term takes 



1 Strictly speaking, Pq(u 1 ,U 2 ) and p (w 1 ,w 2 
constant for non-spherical shapes. 



may not be 



dA 7P -1 



Pi 
Po 



dA XH(— - 
Po 



1 



dA 



fccff 

2 



P_ 

Po 



dA 7p 
Po 



dA 



{2Hf 



dA 



(10) 



where a® = a — 2jpo. The terms J dA 7p + and J dA jp~ 
will be dropped henceforth since we will only consider sit- 
uations where the number of lipid molecules in each mono- 
layer is conserved. 

A couple of important remarks on the physical inter- 
pretation of the above free energy are due here. First, since 
our theory is intended to describe the dynamics of vesicle 
conformation that can be resolved at optical or mesoscopic 
length scales, the free energy should be understood as an 
effective Hamiltonian resulted from an appropriate coarse- 
graining procedure which integrates over fluctuations at 
suboptical scales [^2|^|. The work presented in Refs. [^2| 
p3| makes it conceptually clear that the physical param- 
eters describing the elastic properties of the vesicle mem- 
brane at the mesoscopic length scales, namely, fc e ff, n e fi, 
and A, are "renormalized" parameters, different in princi- 
ple from those parameters describing the elastic proper- 
ties of the membrane at microscopic length scales. For the 
same reason, the constant ag conjugate to the "apparent" 
membrane area under an optical resolution is an effective 
membrane tension p3| , and is distinct both from the mi- 
croscopic surface tension and from the mechanical frame 
tension which is associated with the equilibrium shape and 
which can be accessed by micromechanical measurements 

This issue has not been addressed before in the 
ous works on membrane and vesicle dynamics 
And, one might argue that, when fluctuations are suffi- 
ciently small to allow a Gaussian (linear) theory of the 
fluctuations, the distinction between the effective elastic 
parameters at the mesoscopic length scales and the elas- 
tic parameters at microscopic length scales would only be 
conceptual rather than quantitative, as the conventional 
wisdom on linear theories of fluctuations may lead one 
to believe. However, the work reported in Refs. [f22"| , |23f| 
clearly shows that care must be taken when we apply the 
conventional wisdom to fluid membrane systems. In par- 
ticular, if we consider a vesicle whose membrane is nearly 
incompressible at microscopic length scales and which is 
in either a floppy state or a entropic-tense state the 
effective area compressibility modulus k e s of the vesicle 
membrane at mesoscopic length scales originates almost 
entirely from integrating over small fluctuations at subop- 
tical length scales and is smaller than the microscopic area 
compressibility modulus. We will show in later discussions 
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that our new interpretation of the physical parameters in 
the model free energy bears consequences for the quali- 
tative as well as quantitative characterization of vesicle 
dynamics. 

The second remark concerns the validity of the ex- 
pansion form of the free energy, and relatedly, the choice 
of the reference state, represented by po- As discussed in 
Ref. the effective free energy associated with optical 
length scales involves a nonlinear area elasticity. There- 
fore, the expansion, which only contains linear area elas- 
ticity, is valid in principle only when deviations from the 
chosen reference state po are small. Thus, for our consid- 
erations of near-equilibrium states, a natural choice of po 
is 

Po = N s /A , (11) 

Finally, with definitions tjx^ = p^/po — 1, Eq.([Io|) can 
be rewritten as 



F = dA 



k, 



dA 



«cff 



I 2 



J d A^(2H) 2 + J dA XH(^ 



cr / dA 



(12) 



3 A systematic and general formulation of 
vesicle dynamics 

In this section, we will present a formulation of vesicle 
dynamics. The formulation is both systematic, relying on 
various basic concepts and notions from condensed matter 
physics, and general, in that it can in principle be applied 
to vesicles with arbitrary geometry. 



3.1 Surface hydrodynamics of the membrane: 
Equations of motion 

Following our discussion of vesicle dynamics in Introduc- 
tion, we begin with a hydrodynamic description of the sur- 
face flows within the two fluid monolayers constituting the 
vesicle membrane. The description treats the two mono- 
layers as two individual, but weakly coupled, systems of 
quasi-two-dimensional fluids; it also takes into account the 
fact that each monolayer interacts with its corresponding 
bulk fluid environment. Moreover, it is assumed in the de- 
scription that the motions within the two monolayers take 
place under isothermal conditions, leaving to be consid- 
ered for each monolayer only two intrinsic hydrodynamic 
fields, the density field and the surface velocity field. Con- 
sequently, the surface hydrodynamics of each monolayer 
is principally governed by two equations of motion: one 
states the law of mass (particle) conservation, and the 
other expresses the generalized law of momentum conser- 
vation. 

To facilitate the presentation, we first introduce some 
additional notations related to the dynamics. Let's con- 
sider a fluid surface described by R{u l , u 2 , t), and let the 



trajectory of a particular material particle be represented 
by u a — u a (£, t), a — 1, 2 where £ labels the material par- 
ticle. The covariant components of the intrinsic velocity 
of the particle are then defined by 



_ du a _ du a {£_,t) 
= ~dT = 



dt 



(13) 



The material velocity of the particle in the three dimen- 
sional embedding space can be obtained easily: 



dR(u l ,v?,t) 



dt 

W a t a + d t R 



= W a t a 



dR(u\u 2 ,t) 



Of 



(14) 



U may also be expressed as a decomposition into compo- 
nents tangential and normal to the surface, 



U = U t + U n n, 



(15) 



where U n = U ■ n and U t = U a t a . U a is in general not 
equal to W a as a result of the time- varying geometry, but 
is related to W a as follows, 



U a = W° 



a af} t p ■ (dtR) 



(16) 



For any function f(u ,u ,t) defined on the fluid sur- 
face, a useful time derivative may be defined: 

D t f(u\u 2 ,t)^a-^ 2 d t (a'/ 2 f)=d t f + f^. (17) 



3.1.1 Equations of motion 

We can now write down the four principal equations of 
motion, two for each monolayer, that govern the dynam- 
ics of a fluid lipid-bilayer vesicle. The formalistic deriva- 
tions of the equations are given elsewhere p2] , |2q ]. Given 
the monolayer surface density fields, p ± (u 1 , ir,t), and the 
intrinsic velocity fields, W± (it 1 , u 2 , t), we can define in- 
trinsic particle fluxes, j± = p W±. The subscripts "+" 
and "— " , same as the superscripts, refer to the outer and 
the inner monolayer of the vesicle, respectively. The law of 
mass (particle) conservation for each monolayer can then 
be expressed in covariant form as follows, 



Dtp* + Doi± = 



(18) 



where D t and D a have been defined previously. 

The generalized law of momentum conservation can 
also be written in a similar way. For each monolayer, 
a three-dimensional vector representing the surface den- 
sity of monolayer momentum can be defined by J± = 
mp ± U ± = mp ± (W±t a + dtR), where m is the molecular 
mass of the constituting lipid. The equations of motion 
governing the changes of J± take the following form: 



D t J± + D a (jlJ ± ) = f 



(19) 



where j±J± are the covariant momentum fluxes. / rep- 
resent the total forces per unit area acting on the mono- 
layers. 
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3.2 Forces acting on monolayer surfaces 

In principle should contain several forces of different 
origins: any dissipative forces associated with the intrinsic 
surface shear and dilational viscosities of the monolayers; 
ft, forces which the two monolayers exert on each other; 
ft, the mechanical (restoring) forces associated with the 
free energy Eq.(|l^); and, finally, T , the hydrodynamic 
forces exerted on the monolayers by their corresponding 
bulk fluids. In the following, we discuss each of the four 
types of forces in turn. 

In our theory of vesicle dynamics, we have neglected 
the intrinsic dissipative forces within each individual mono- 
layer, based on the evidence and arguments given in Ref. 



16|. The intermonolayer forces, may be modelled as 



ft = Tb(U 



u: 



ft, 



(20) 



where 



ft = Tb(U+ U~ t 



(21) 

phenomenologically describes dissipation associated with 
the relative motion between the monolayers. Thus, b may 
be called the "intermonolayer friction coefficient." This 
description was originally proposed by Evans et al. 
ft, with = — /~ , represent the normal components of 
the intermonolayer forces, which constrain the motions of 
the two monolayers in the normal direction. 



3.2.1 Mechanical restoring force 

In deriving from the free energy given in Eq.([l2]) the me- 
chanical restoring force acting on each individual mono- 
layer, we may imagine the following infinitesimal virtual 
variation in the geometry of each monolayer: 

R{u\ u 2 , t) -> R{u\ u 2 , t) + 5R{u\u 2 ) , (22) 

and define the restoring forces ft through 

SF ± = jdA (-ff s ) ■ 6R . (23) 

where SF + and 6F~ are the variations in the free energies 
of the two monolayers induced by the shape variation 
6R. 

Two issues are implicit in Eq.(p3[), which need to be 
addressed clearly. The first is the issue of how to evaluate 
F ± . This issue has not been discussed at all in the current 
literature on dynamics of lipid-bilayer membranes. Our 
proposal is to divide F into two "monolayer parts" , F — 
F + + F~ , where 

F ± = j dA ^ ■ (</,±) 2 ± I dA XHcf> ± 

+ l -{jdA^f(2H) 2 + a jdA^ (24) 

is associated with a single monolayer, and then to derive 
the corresponding restoring force from the variation of the 



monolayer free energy. This definition is consistent with 
the bilayer-composite structure of the vesicle membrane. 

The second is the issue of how to deal with the vari- 
ations in the monolayer density fields p + (u 1 ,u 2 ,t) and 
p~ (u 1 , it 2 , t) in their relations to the shape variations of 
the monolayer surfaces. Our proposal ^ concerning this is- 
sue is that the number of molecules associated with any 
local area element dA of a monolayer should be conserved 
under the shape variation of the monolayer. In other words, 



5{dA p d 







(25) 



The reason for this lies in the basic thermodynamic con- 
sideration that the variation in the free energy should be 
equal to the mechanical work done against the restoring 
forces. Briefly, in considering the hydrodynamics of the 
membrane surface, we assume that the thermodynamics 
of a local monolayer element with area dA and num- 
ber of molecules is described by the local free energy 



F A [T.X\.dA[R).R)=dA \ 



h± ) 2 + a i 



. dA {^f(2H) 2 ±XH ( f> ± } 



(26) 



Clearly, the change in the local free energy, which can 
be identified with mechanical work, must be the change 
under constant T and constant N A . 

The mathematical derivation of the restoring forces for 
a vesicle of arbitrary shape is similar to that described by 
Jenkins [^7| and is somewhat lengthy. We will only state 
the final result here: 



ff s =1 - K e s [H (2H 2 - 2K) + AH] - n ± (2H) 



+ {2H 2 - K)^ + 2H 2 



n 



(27) 



- V t ^TA(l 



where A = (1 / \fa)d a {a a ^ y/adp) is the Laplace-Bcltrami 
operator, V t = t^a a ^d a is the gradient operator defined 

on the surface, and 7r± = — — + fc c ff/ ,± H — ^j— ( < / ,± ) 2 may 

be considered as monolayer surface pressures under the 
condition of planar geometry. 



3.2.2 Hydrodynamic forces due to the bulk fluids 

The hydrodynamic forces exerted on the monolayers by 
their corresponding bulk fluids can be determined by solv- 
ing the equations of motion for the two bulk fluids sepa- 
rated by the membrane under appropriate boundary con- 
ditions. 

Equations of motion 

Since the bulk solvents for lipid-bilayer vesicles are always 
aqueous, it is a good approximation to consider the flu- 
ids as being incompressible. Let the velocity and pressure 

2 Our proposal appears to be the same as that used in 
Ref. H. 



Ling Miao et al.: Dynamics of shape fluctuations of quasi-spherical vesicles revisited 



7 



fields in the two bulk fluids be represented by v^°\r,t), 
v^(r,t), and p(°'(r,t), p^(r,t), respectively, where the 
superscript "o" refers to the bulk fluid outside the vesicle 
and in contact with the (+)-monolayer and the superscript 
"i" refers to the bulk fluid inside the vesicle and in con- 
tact with the (—) -monolayer, ft is then straightforward to 
write down the hydrodynamic equations: 



V • v^l = , 



(28) 



By quasi-spherical vesicles, we refer to those vesicles with 
very small excess areas such that their equilibrium shapes 
are spherical. Thus, both equilibrium shape fluctuations 
and non-equilibrium shape deformations may be consid- 
ered as small perturbations around the spherical equilib- 
rium shapes. Under such conditions, an analytical render- 
ing of the theory becomes possible, provided that certain 
approximations are made. 



pb 



T + f>« ■ V)v^\ = Vp( a > + rjV 2 v^ , (29) 4 j Q uasi . spheri cal vesicles: specific notations 



where pb and r\ are the mass density and the shear vis- 
cosity of the bulk fluids, respectively, and the superscript 
a = i, o. 

Boundary conditions: kinematic matching 

Each monolayer of the vesicle provides a boundary surface 
for the corresponding bulk fluid. Boundary conditions, 
which specify V (a) (u 1 , u 2 , t) = v^(r = R^u 1 , u 2 , t) , t) , 
are required for completely determining solutions to the 
above equations. Under very general circumstances, V^ a ' 
may be different from the flow velocity in the correspond- 
ing monolayer, J7 ± (w 1 , u 2 ,t). In the absence of any practi- 
cal evidence for the general scenario, we choose to match 
the bulk kinematics at the boundary surfaces with the 
monolayer surface kinematics, i.e., 

V (o) (v\tt a ,t) =U + (u\u 2 ,t) , (30) 

^(ttWl^r^ 1 ,)!^). (3i) 

Bulk hydrodynamic forces 

Once the complete solutions to the bulk hydrodynamic 
equations are obtained, the corresponding hydrodynamic 
stress tensors, TW, are readily derived as, 

T^(r,t) = _pW( r ,t)I+^ \w^(r,t) + (W a) (r, t)) T ] , 

(32) 

where I is the unit tensor, and the superscript " T " indi- 
cates the transpose of a tensor. The forces per unit area 
exerted on the two monolayers by the corresponding bulk 
fluids can finally be evaluated as follows, 



T+(u\u 2 ,t) = T (o) {r,t) 
T-(u\u 2 ,t) = - T«(r,i) 



r=R(u\u 2 ,t) 
r=R(u 1 ,u 2 ,t) 



n(u 1 , u 2 , t) , 



■ n(v},u , t) 



where n(u 1 , u 2 , t) indicates the local normal vector of the 
membrane surface that directs towards the exterior of the 
vesicle. 



4 Dynamics of shape fluctuations of 
quasi-spherical vesicles 

Having formulated in the previous sections a systematic 
and general theory for near-equilibrium dynamics of lipid- 
bilayer vesicles, we will in this section present the appli- 
cation of the theory to systems of quasi-spherical vesicles. 



For describing the surface geometry of a quasi-spherical 
vesicle in the embedding three-dimensional space, the most 
convenient coordinate system is the spherical coordinate 
system (9, (p, r), with local unit basis vectors e#, e v , and 
e r . The quasi-spherical geometry of the membrane surface 
of the vesicle can be represented by 

R(0, (p, t) = R(9, <p, t)e r = Ro[ 1 + u(9, ip, t) ]e r , (33) 

where Ro is the radius of the spherical equilibrium shape, 
and u(9, ip, t) describes an arbitrary, small shape perturba- 
tion. The two local tangent vectors on the quasi-spherical 
surface are then given by 

_ dR 
o9 



dR 

t v = -Q— = R sin 9e v + R v e r , 



(34) 



where Rg = dR/89 = Ro du/89 and R v = OR/ dip = 
Ro du/dip. 

It follows from Eq.([l6|) that for this geometry, the co- 
variant components of the velocity of each monolayer can 
be found as 

E/| = W± + a al3 {t ■ e r )d t R = W| + a al3 R p d t R . 

Clearly, in the scheme of linearization where quantities are 
expressed accurate only to first order in the perturbation 
represented by u, 



Uf = U£t a = Rq W e ±e e + (R sind) 



(35) 



The components of U± — U^eg + U^e ip + U^e r will also 
be needed later and are expressed here as well to first order 
in the perturbation 

ut(9,<p,t)=RoW£(e,ip,t) , 

U±(9,p>,t)=R sm9W£(9,p,t) , 
U±(9,<p,t)=R u(9,ip,t) , (36) 

where we have introduced a short-hand notation, g(t) — 
dg(t)/dt, for the time derivative of any function g(t). An- 
other quantity to collect here is D a W±. Within the scheme 
of linearization, 



fdWl 



■ cotewi 



dip 



(37) 
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Two differential operators will be defined here for later 



use, 



(38) 



rr d Id 

ad sin dtp 

and L? = — Vl ■ Vl- Clearly, by invoking I? on the spher- 
ical harmonics, 3^ m , we have 



L 2 y tm = t[i + i)y tr , 



(39) 



Finally, we write down the expansions of all the rel- 
evant fields, u(9,(p,t) and ± (0, tp,t) in the basis of the 
spherical harmonics, yi m (9, ip), 

u(o, <p,t) = y j u lm (t) y em 

m=l 
1=2 m=-i 

m=t 

= <+EE 1>L(t)y im , (4i) 

1=2 m=-i 

<j) S (8, if, t) = <$ + Y,*fJt) 

£,m 

rn=£ 



1=2 m=- 



where two alternative fields 



4> A (9,<p,t) 



_ <l>+(0,<p,t)-<l>-(fi,(p,t) 



<P(0,<p,t) = 



_ 0+(^vJ,i) + 0"(».V,*) 



(43) 
(44) 



are used instead of 4> ± (9, ip, t) for later convenience. <j> A = 
(0o — 0o~ an d <j)§ = (0q + <j> ~)/2 — are constants 
characterizing the equilibrium state of the vesicle; and 
{v>e m (t), ipfrn(t), tpf m (t)} represent the perturbations. 

The reasons for not including I = 0, 1 modes in the 
above expansions may not all appear obvious. The I = 
mode is excluded in Eq.(|40|) because we will not be con- 
sidering this mode. This is related to the fact that the 
dynamics of this mode in the fluctuation spectrum is not 
only slaved by the dynamics of the other modes, but is also 
unresolvable experimentally. The t = 1 mode in the per- 
turbation expansion corresponds to a simple translation 
of the vesicle and its amplitude is thus set to zero with- 
out loss of generalit y. T he exclusion of the t = 1 modes 
in Eq.([4l"l) and Eq.([42|) is due to the fact that they are 
not coupled to the shape changes of the vesicle within the 
scheme of small perturbations. 



4.2 Approximations 

In the application that will be presented in the follow- 
ing, the most significant approximation which has been 



made is that the near-equilibrium dynamics of the vesicles 
are purely dissipative, or overdamped. In other words, the 
inertial effects in both the bulk hydrodynamics and the 
monolayer surface hydrodynamics are neglected. This ap- 
proximation has apparently been made in all theoretical 
work on vesicle dynamics, based on the estimates put for- 
ward first by Milner and Safran . Underlying the consid- 
eration of Milner and Safran was, however, the assumption 
that the only relevant dynamics of a membrane vesicle was 
the dynamics of the membrane conformation. Therefore, 
although we adopt the same approximation in this paper, 
we feel that additional considerations to those that led 
to the earlier estimates are needed, given the fact that in 
addition to the surface geometry field monolayer density 
fields are also being considered. Thus, we briefly state this 
issue here. 

Let's consider the bulk hydrodynamics first. Making 
the approximation amounts to neglecting both the iner- 
tial term and the non- linear convcctivc term in the Navier- 



Stokes equation written in Eq.(28). By estimating the re 



spective ratios of the two terms to the viscous-force term, 
we may have a guideline on whether the approximation 
is justified or not. The ratio of the inertial term to the 
viscous-force term can be measured roughly by a dimen- 
sionless quantity 



K 



p h L 2 p h \dv/dt\ 



r)to 



7]\W 2 V\ 



(45) 



where L and £o are j respectively, the typical length and 
time scales characterizing the spatial and temporal varia- 
tions of v. Another dimensionless quantity measures the 
ratio of the convective term to the viscous-force term, 

= p h vL ^ p b |(t) • V)v\ 



V 



T]\V 2 V\ 



This quantity is the conventionally defined Reynolds num- 
ber pq |. 

In the case of vesicle dynamics, motions in the bulk 
fluids are induced by motions of the vesicle membranes, 
which may be considered as a linear composition of modes 
of motion characterized by different wavelengths. For each 
given mode L = a p,m should correspond to the charac- 
teristic wavelength of the mode. For a vesicle of 20/im 
radius, a may range from 1 to 60, covering those modes 
that can be resolved under optical microscope. For a self- 
consistency check of the no-inertia approximation, the time 
scale to should represent the shortest time scale character- 
izing the overdamped dynamics of a vesicle, rather than 
simply the time scale associated with the overdamped re- 
laxation of a pure bending mode, as was the case con- 
sidered by Milner and Safran jl]]. For pedagogical reasons, 
we will leave the quantitative self-consistency check to the 
discussion, where more understanding of the various relax- 
ation time scales will have become available. 

To estimate the order of magnitude of 1Z, we replace v 
by I /to, where I should be a measure of the shape deforma- 
tions of the vesicles from their equilibrium shapes, and is 
always much smaller than L in the cases that we are inter- 
ested in. Consequently, we have TZ/TZ = l/L <C 1. Based 
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on this analysis, we may always neglect the convective 
term in our considerations of near-equilibrium dynamics 
of vesicles, provided that 1Z is smaller than 1. Equations 
( fj8| ) thus become 



V • v (a) = , 
Vp (a) = ?7VV a) 



(46) 



Similar order-of-magnitude analysis can be made re- 
garding the approximation of neglecting the inertial and 
convective terms in the two-dimensional counterpart of 
the bulk Navier-Stokes equation, Eq.(^9|). The order of 
magnitude of the surface inertial term is given by p m \U\/to, 
where p m is the surface mass density of a lipid monolayer, 
and the order of magnitude of the forces acting on a mono- 
layer may be represented by |T| ~ r\vjL ~ rj\U\/L. The 
ratio of these two terms is then given by 



1l m = 



K . 



Ph 



(47) 



It is easy to work out that lZ m <C 1Z. Thus, The ap- 
proximation that the monolayer surface hydrodynamics 
is overdamped is well justified. Under this approximation, 
Eq. (|l9|) reduces to 



/± = ft +T ± Tb(U 



u: 



o 



(48) 



4.3 Mechanical restoring forces 

The general expressions for the mechanical restoring forces 
given in Eq.(B7ft should be evaluated now for the specific 
case of the quasi-spherical vesicle described by Eq.(33) to 
first order in the perturbation fields. 

In the scheme of linearization, the mean curvature H 
and the Gaussian curvature K take on the following forms, 

U ZHo c — 



I .11, 



K= ~k + ~k^ {£ l)Ulmyim 



(54) 



The Laplace-Beltrami operator A and the surface gradient 
operator Vt become 



i 

Ro 



if <P(8, tp) is any function of first order in the perturbation 
fields. 

By the use of the above simplified expressions, the sum 
of the two monolayer normal components is evaluated to 
be 



It is more useful to decompose the above vector equa- 
tions into their components in the normal and the tangen- 
tial directions. Moreover, it turns out that, in the normal 
direction, only the sum of the two monolayer equations is 
relevant, which reads 



(/+ + /") -n+(T + +T-)-n = Q 



(49) 



(fls + /„) ' n 
2? 



4fc off 



2t >-V r K ctf <±rv eS E 

R:-l^\-m Etu ^-—^ m 



i.m 



Rq 



4^ + ^(^ + 2)^-1) 
Ho K n 



X 



(55) 



In the tangential directions, the more convenient expres- 
sions can be obtained by the difference and the sum of the 
two monolayer equations, which look like 

(/rt.t - /rs.t) + (Tt - T~) - 2b (U+ - C/ t ") = , (50) 



J IS 



(T 







(51) 



In principle, the last three equations are applicable to 
vesicles with arbitrary shapes. In applying these equations 
to the case of a quasi-spherical vesicle, we will only con- 
sider contributions which are first order in the perturba- 
tions. The first-order contributions from both ff a and 
can be obtained exactly without approximation. The situ- 
ation concerning is not as well controlled, since an ex- 
act solution would involve solving the bulk hydrodynamic 
equations for boundaries of arbitrary shapes. A practi- 
cal approach, which is consistent with the linearization 
scheme, is to make an approximation, 

T ± =T±;, (52) 

where Tq is obtained by solving the bulk hydrodynamic 



equations given in Eq.(46) for spherical boundaries that 
coincide with the spherical equilibrium shape of the vesi- 
cle. 



where 



and 



with 



Co 



Rq 



E t = (£ + 2)(£-l) 



£(£+l)+r ^- 



TO 



Mi 

Rq 



(56) 



(57) 



(58) 



The two tangential forces, t and / rs t , are, to first 
order in the perturbations, given by 



/rs,t 



-fceff(l + ^)V L ^ T A(l + ^) V L ff 



The really relevant quantities are not the tangential forces 
themselves, but the following ones, 



Vl • (/+t 
A 



i?3 



/™,t) = E^ +1 ) 

l,m 

(£ + 2){£-l)u im + 



2^55" Vim 



Rl 



(59) 



y<: 



10 



Ling Miao et al.: Dynamics of shape fluctuations of quasi-spherical vesicles revisited 



and 



VL-(/+ t +/- t )=E^ +1 



Z r>2 Vim 



n 2 



4. 



( I/I 



(60) 



where and p$ are the hydrostatic pressures in the 
bulk fluids outside and inside the vesicle, respectively. 

When the tangential components, T^ t are concerned, 

it is more convenient to use Vl ■ T$ t , which look like 



Vi, • T2 



4.4 Forces due to the bulk hydrodynamics 

The forces, T^, defined in Eq. (|52|) can be obtained based 
on an adaptation to our problem of the classical Lamb 
solution [g9| . This approach has already been used in the 
previous work on vesicle dynamics by Schneider et al. p0| , 
Yeung and Evans jl6|], and Seifert |Q. We will, therefore, 
only summarize the final results here. 



R 



1=2 m=-l 



R 



L E E [(^2)^ + (2f+l)^ 

'0 fi O fi 

Vl • T " t = 

-i)x£ + (2£ + i)y« 



y 



lm 



u i> n /> 



2 m=-( 



(66) 



4.4.1 The Lamb solution 



4.4.2 Surface continuity equations and kinematic matching 



In the Lamb solution, an alternative form of the boundary 
conditions, 

t> (a) (r, 9, (f,t)\ r=Ro = V {li \d,ip,t), is adopted, where the 
following three scalar quantities are used instead of the 
three components of (9, ip, t): 

X^(9,^,t) = V {a) (9,^t)-e r , 

y W (9, ip, t) = -Ro V • V (a) {9, <p, t) , (61) 

Z^{6,ip,t) = R e r • (V x V (a) (9,(p,t)) . 

The above three quantities can be expanded in the basis 
of the spherical harmonics, 



x^\9,^t) = Y J x^(t)y im , 



lm , 



(62) 
(63) 
(64) 



The normal and the tangential components of T Q can now 

be expressed in terms of X$*(t), Y${t), and Z^(t). The 
normal components read as 



T, 



0,n 



-Po 



Rq 



V 



=2 m=-l 



(i) V 



^O.n — Po 



u 1=2 m=-£ 

- V 

Ro 



1 



3 V 



EE 



lm i 



(65) 



The "amplitudes" X&>, X®, Y^ and Y® can be di- 
rectly related to the amplitudes of the three perturbation 
fields, u im (t) and ^f m (t) as defined in Eqs.©-@. The 
connections are provided by the two surface continuity 
equations given in Eq. ( |l8|) and by the kinematic-matching 
conditions given in Eq.(|30|). Expressed to first order in the 
perturbations, Eq.([l8) read as 



Po 



1 



^ (<?,¥>,*) -2u(6,<p,i) . (67) 



1±, 



Substitu ting the kinematic-matching conditions Eq.(30) 
into Eq.(|u2|) and Eq.(|6^), using the expressions for 
given in Eq. (|36|) , we arrive at 

X W (9, ip, t) = X® (9, if, t) = Ro u(9, <p, t) , (68) 
Y {o) {9 1 ip,t) = -2Rq u(9,ip,t) (69) 



— Ro 



Y®(9,tp,t) = -2R u(6,tp,t) 



(70) 



Rq 



1 d 
^9 89 



(W e _ 



sm f 



8W_ 

dip 



Applying Eq. (^) to the last two equations, we finally have 



X^(t)=XZ(t) = Rou em (t), 



r(0),.S RO 



1=2 m=- 



Y^(t) 



Ro 



~1 i>tnif) . 
°0 



(71) 
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By substituting Eq.fln]) into Eq.@ and Eq.@, we tution of both Eq.(||) and Eq.(j72|) into Eq.Q yields | 
arrive at the following relevant expressions: 



<P { o ] 



V 



£(£ + !) 



(4£ 3 + 6£ 2 -Al- 3)% 



3[l + (2l+l)^; 



3(2£+l + <^) )E 



(72) 



V L .(T+ t -T - t ) = 



2(21+1)0* ]s ■ 
7 L .(T + t + T - t ) 



Ro 

ytm 

= A 

Rq 



(73) 



2(2£+l)^ - A 2(2£+l) ]E 



(2£+l)ut m 

yim ■ 



(74) 



£{£+!) 



(4£ 3 + 6£ 2 -4£- 3)u tr< 



3[l + (2* + l)^] • A , 3(2^+ 1 + 0*) )s 



-Ui-n 



1 

KcffE j 

"If 

4fc c ff ,£ 



tf) 2 



4^ + 4(/ + 2)(/-l) 

-Kq -Kg 



(77) 



Similarly, combining Eq.@, Eq.@ and Eq.Q based 
on Eq.(|5C|) leads to 



£(£ + !) 



2(2£ + 1) 



4&i? 



-^{£ + 2)(£-l)u em -2^i> 



Rq 



2k eS ^ A 



Rq 



Vim i 



(78) 



and combining Eq. (|6C|) and Eq.(|74j) according to Eq.(|51 
gives 



4.5 The intermonolayer frictional forces 



The evaluation of the intermonolayer frictional forces ap- 
pearing in Eq.(pO 



is straightforward. By applying Vl on 
Uf 1 given in Eq.(35|), we obtain, to first order in the per- 
turbations, 



V L • Uf = R 



f dwj 



Ro {D a W c ± 



cot ewi 



Ro 



dWl 

dip 



1 ± 



2ii 



(75) 



where Eq.([37|) and Eq. (p7|) have been used. It follows im- 
mediately that 



) (76) 



4.6 Equations of motion 



A set of 3 linear equations of motion can now be obtained 
for the three dynamic variables characterizing a particular 
mode of perturbation, ue m (t), ^^(t) and tp£ m {t). Substi- 



£{£ + !) 
2{2£ 



{21 + l)u im - 
-1) 



2(2*+!)^^ 



1 - 



^ 2 



fin 



Mi 



R 



r (f+2)(l-lK m -^t 

Hq 



Hq 



(79) 



In principle three independent modes of dynamics can now 
be determined by solving the three equations of motion 
given above. 

A closer analysis of the equations shows, however, that 
only two modes of the three are relevant on the time scales 
that are experimentally accessible, due to the fact that 
there is an inherent separation of time scales involved in 
the problem. Three basic time scales may be defined based 
on Eq. (0) , Eq. @ , and Eq. © , 



tjR% £ 3 
K eff £(£+l)E e 
V R {2£+l) 



tA 



2bR 2 



cir 



£(£ + !) 



k cS £(£ + !) ' 

An order-of-magnitude estimate based on 

erg 



10 



b = 5 x 10 



12 erg , 
7 erg • s 



30 



i?o = 10 /J.m , 



(80) 



(81) 



rpRl 



= 10 



3 The equation at the zeroth order is given by p^ — p^ 



2f/R . 
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yields for £ = 2, 

t c ~0.2s, t A ~0.5s, ts-2xl0 _7 s 

Clearly, the relaxation of ipf m , which is characterized by 
tg, is much faster thi 
Thus, we can assume 



ts, is much faster than the relaxations of ui m and tphn- 



4>Fm * , (82) 

over the time scales that characterize the relaxations of the 
slower modes. Using then Eq.([79]) to eliminate ij)f m from 
the other two equations, we arrive finally at two relevant 
equations of motion, written in the following in a matrix 
form 



7m 



An v oB 



2x2 D 2x2 



Ufa 
'fr 



= -c 



2x2 



Ut: 
i m 



(83) 



The two coefficient matrices are defined as 



1 



and 

B2x2 = 



V 

1(1+1 
to 



-3 + 0t4tf - 



46 + 2ft(l + (^ 
l-(0if) 2 



A\2 



\ 
J 



(84) 



-E t + 2<f£\t 
-A, [l-Wtf) 2 ] 



2k [1 



(85) 



where to = t/Rq/ k c s defines a time scale, and 



bRo 



XR 

Keff 



k = 



fc c ffi?Q 
Keff 



are dimensionless parameters, and where a few short-hand 
notations have also been defined 



r e = {21 + 1){2£ 2 + 2£- 1) 
\t = Xo(£ — 1)(£ + 2) . 



fa = 21+1 



(86) 



Further approximations can be made in the analysis 
of Eq. (|8"3|) . The dimensionless parameters bo and ko ac- 
quire rather large values for those values of the physical 
parameters quoted in Eq.(j3l|), specifically, b Q w 5 x 10 6 , 
and ko ~ 3 x 10 7 , whereas parameter <p^ is expected to 
be much smaller than one. Thus, both A2 X 2 and B2 X 2 can 
be simplified as 



Atx2 — 







-4b 



Bo 



£{£ + !) 



tr 



-Et 
-Xe 



24>qXi 



Ac 
2k 



(87) 



This yields a rather simple form for matrix C2X2, 
t(t + 1) ( Ab (Ei - 2(f>£X e ) -4b Q X ( 



C?x2 — 



4b r e t 
_ ( Cn C12 
C21 C22 



-A^-Tj 



2k a r e 



(88) 



A comment is worth making here, which should help 
clarify the connections between our theory and the pre- 
vious theories on the vesicle dynamics |Q,^,|l6|. It con- 
cerns the two approximations made by use of Eq. J82] ) and 
Eq. (^) . The first approximation is actually equivalent to 
imposing the so-called "local area incompressibility con- 
straint" that has been used in some of the previous works 
The physical reason underlying the approxima- 
tion, expressed in Eq. (|82]) , thus provides a rationale for 
the use of the constraint. However, it is incorrect, in our 
opinion, to interpret in general, as the authors of Ref. [ [l6| 
did, this constraint as representing incompressibility of the 
monolayers as fluids, to be understood in a way similar 
to the concept of incompressibility of bulk fluids, since 
the constraint eliminates only one of the two mechanisms 
by which the monolayer density fields can change - the 
mechanism through changes of local surface geometry of 
the monolayers. The second approximation implies that, 
as far as the induced hydrodynamic motions in the bulk 
fluids are concerned, the surface flow effects that are asso- 
ciated with the overall compressibility of the monolayers 
are negligible. This approximation has the same effect as 
the approximation used in Ref. fl6| (the second equation 
in Eq. (A. 6) therein). 



4.7 Dissipative dynamics of shape fluctuations 

Finally, two dispersion relations, which characterize the 
two independent dissipative modes of the vesicle dynam- 
ics, can be obtained as the two eigenvalues of C2X2, 

1 



tt±(£) = -[f2 c {£) + Q A {£)} 



1± 4/1 



f2 A (£) ■ qjj) 
\o c (£) + n A (£)} 



(89) 



where 

^W = c 22 

O c (£) = C n 



£(£ + l)fceff 



2bRl 



£{£ + !){£ -!){£ + 2) 



vR 3 r e 

£(£+!)(£ -!){£ + 2) 



[£(£ + 1)k cS + E a 
X 2 

£(£ + 1)/Sa + E a + — 
k c s 



In the above expressions, the definition 



E„=Rl[ao-k e -< 



has been used, and 



2k 



c rr 



(90) 
(91) 

(92) 



defines a new bending rigidity, which will be called "ap- 
parent" bending rigidity. 
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At this stage, the quantity of our main concern, namely, 
the time-dependent correlation function of shape fluctua- 
tions of a quasi-spherical vesicle, (w| TO (t)u^ m (0)), may be 
written generally as 



\Utr, 



(93) 



Clearly, /3_ + /3+ = f . To obtain analytical expressions 
for the amplitudes, /3 T , we employ Onsager's regression 
hypothesis that the equations of motion governing the 
macroscopic near-equilibrium dynamics can be applied to 
the dynamics of spontaneous fluctuations Following 
the formalism given in Ref. we thus have 

(u* em (t)u em (0)) = (kml 2 ) (e _C *)n+(t*J m ^> (e~ Ct h 



(94)' 



which leads to the following specific expressions for (3^: 



f2 + 




J7_ 




i 












i 








Q- 




1 




n+ 




i?_ 



Qa — — C21 
Qa - (2- + f2 A 
J? + — Qa + C21 
i7+ — [2 a — [2 a 



(\u£ m \ 2 ) 

X(£ -!)(£ + 2) 



2fc c ffi? 

{\ut m \ 2 ) 

A(l- l)(l + 2) 
2k c gR 



(95) 



(\ u em\ 2 ) an( i k^im^tm) represent some of the static cor- 
relation functions, the calculations of which are briefly 
sketched in Appendix A. 



5 Discussion 

The dispersion relations, Eq.([S9"|), predicted by our the- 
ory appear formulistically identical to those worked out 
by Yeung and Evans jl6| and by Bivas et al. jlj. But, 
our definitions of two of the three quantities involved in 
the dispersion relations, namely, f2 c (£) and f2 & (£), differ 
in detail from those previous results. More importantly, 
we feel that interpretations of those physical parameters 
contained in the dispersion relations, K c ff, k c g , A and E a , 
need to be reconsidered. 

We have already remarked on the issue of interpre- 
tation in Section 2.2, when introducing the effective free 
energy describing a vesicle with a coarse-grained configu- 
ration characterized by i^u 1 , u 2 , i), and j o ± (w 1 , it 2 , t). As 
we have not yet performed for systems of quasi-spherical 
vesicles the coarse-graining procedure similar to those de- 
scribed in Refs. p2|p3[ , we are only able to give a semi- 
quantitative discussion of the effective mesoscopic param- 
eters. 

If we consider systems of quasi-spherical vesicles where 
fluctuations about equilibrium shapes are small, where 
a Gaussian theory would suffice in describing the small- 
scale fluctuations, we may argue that the renormalization 



of bending rigidity due to nonlinearities ]32| can be ne- 
glected and that K ff ~ K m — A 2 n /2fc m , where n m , A m , 
and fc m represent the phcnomenological parameters corre- 
sponding the microscopic cut-off length scale. Note that 
the part subtracted from K m arises from the coupling be- 
tween membrane curvature field and the density-difference 
field But, making a quantitative statement about K ff 
is not trivial, at least, not straightforward, even from the 
point of view of experiments. In the so-called flicker-noise 
analysis experiments, a bending rigidity n and a "surface- 
tension" constant S e are inferred from fitting the experi- 
mental data on (|tt^ m | 2 ) onto a functional expression, 



fc R T 



{£-!){£ + 2) K [e{£+1) + S c 



(96) 



Based on our theory and calculations presented in Ap- 
pendix A, which yields 



\Utn 



k B T 



(t-l)(l+2) Ka [t(l +!) + £} ' 



where 



S : 



(ToRl - /c c ffi? 2 (</>jf ) 2 + 



il 1 

k c tf 



(97) 



(98) 



the experimentally obtained bending rigidity k should be 
identified not with k c k , but with K a , the apparent or renor- 
malized bending rigidity. 

k e $ is expected to be lower than its microscopic coun- 
terpart due to the renormalization effect when the vesicles 
are either in or close to the floppy states, as we have al- 
ready mentioned in Section 2.2. It is our opinion that the 
values of area compressibility moduli for various single- 
component membrane systems reported so far in the lit- 
erature |^,^3| are actually the values of the microscopic 
area compressibility moduli, as these values have been 
obtained when the membranes under observation are in 
stretched-tense states [23] . In the absence of any available 
experimental data on k e g, we make an estimate based on 
the semi-quantitative derivation of the parameter given 
in Ref. |23] by using a short-distance cut-off of optical 
lengths, Ji 1 = 100 nm and obtain fc e ff — 30 erg/cm 2 . 
This value is lower by 2-3 fold than the values of the mi- 
croscopic area compressibility moduli typically quoted in 
the literature . 

The parameter, A, is a rather elusive one, due to the 
fact that understanding of the physical origin of A m is 
limited. We will follow the discussion on it put forward 
by Seifert and Langer [l5]| and put an estimate on A m as 
A m = 2k Tn d, where d ~ lnm, representing roughly half 
of the monolayer thickness. Taking d ~ 1 nm and fc m = 
100 erg/cm 2 , we have A m ~ 2 x 10~ 5 erg/cm. It is not 
difficult to see, if one follows the kind of coarse-graining 
procedures described in (22), that A ~ A m . 

The parameter, E a , should according to our calcula- 
tion be related to the experimentally inferred parameter 

A 2 

E a = n a S c - — , (99) 
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as, based on Eq. (|96|) and Eq.(|97|), S should be identified 
with S c . 

Our analytical expressions for (3^ differ also qualita- 
tively both from those given in Ref. Jl6| and from those 
given in Ref. jl7| . Explicitly, the qualitative corrections to 
the expressions given in Ref. [[l6| are 

A _ ± (K»| 2 ) n A / A(i-i)(i + 2) \ 2 
T o+-0- \ 2k cS R ) 

Quantitatively, however, the corrections turn out to be 
negligible, as a quick, order-of-magnitude estimate by us- 
ing the values of A and k c R quoted above and Rq = 20 fim 
shows. But, the corrections to the expressions obtained in 
Ref. jlTj become significant quantitatively, as we will show 
numerically in the following. 

Based on the above discussion, we now present a sys- 
tematic numerical analysis of the analytical results given 
in the previous section. The analysis itself is obviously 
straightforward, but it brings out some numerical conse- 
quences that are very relevant to experimental studies of 
the dynamics of vesicle shape fluctuations. We first ana- 
lyze two cases, where two important parameters take on 
different numerical values: case a) k e g — 100 erg/cm 2 , 
b = 10 7 erg • sec/cm 4 , and case b) k c s = 30 erg/cm 2 , and 
b = 2 x 10 8 erg • sec/cm 4 , whereas the other parameters 
assume the same values in both cases (see the captions for 
the tables and figures). Table 1 together with Figure 1, and 
Table 2 together with Figure 2 illustrate the numerical re- 
sults, corresponding to case a) and case b) , respectively. 

As the numerical results clearly show, case a) and case 
b) yield two rather different scenarios for the relaxation 
dynamics of shape fluctuations of quasi-spherical vesicles. 
In case a), where the values of fc e ff and b are chosen to 
conform to the values quoted canonically in the literature, 
the time scales characterizing the two modes are separated 
by almost two orders of magnitude, and the slower mode 
takes up almost all of the full amplitude. Given the typ- 
ical temporal resolution of milliseconds and the typical 
accuracy of 10% in determining fluctuation amplitudes in 
standard flicker- noise analysis experiments j34|, only the 
slower mode can be effectively resolved by experimental 
observations and analysis. The different physical nature 
of the two modes are also illustrated by Figure 1. The 
dispersion relation describing the slower mode is well ap- 
proximated by f2 a (l), whereas that describing the faster 
mode is well approximated by i7^(/). The functional form 
of J7 a (/) is the same as that of the dispersion relation de- 
rived by Milner and Safran for the relaxation of a pure 
bending mode in the absence of the coupling between 
membrane geometry and monolayer density fields j?J . The 
consequence of the coupling is, however, present in our 
result. It is indicated by the fact that the bending rigid- 
ity appearing in f2 a (l) is the apparent one, K a , which is 
the result of renormalization of K e ff by fluctuations in the 
monolayer density fields in the presence of the coupling. 
The reason for the renormalization effect is also revealed 
in Fig. 1: Fluctuations in the difference of the monolayer 
density fields relax much more quickly than pure shape 
fluctuations, i.e. £2a{1) ^ ^ c (0- 



10000 




Fig. 1. Various dispersion relations relevant for describing the 
two dissipative modes of vesicle dynamics for K a = 1CF 12 erg = 
24k B T, r\ = 0.01 ergs/ cm 3 , A = 2 x 10" 5 erg/ cm, E = 10, 
R = 20 /im, fc fj = 100 erg/ cm 2 , and b = 10 7 ergs/ cm 4 . 



In case b), the quantitative differences between the two 
relaxation time scales are smaller than in case a) , and the 
slower mode takes up a much smaller fraction of the full 
amplitude. Relevant to experimental studies is that, for 
the first few low £ values, both relaxation modes can be re- 
solved with the typical experimental resolutions. Another 
point to note is that there is a reversal in the physical 
nature of the two modes in comparison with case a), as 
illustrated by Figure 2. It is now the faster mode that de- 
scribes approximately the relaxation of bending deforma- 
tion, where the governing bending rigidity is, however, K e s 
instead of K a . In other words, the renormalization effect is 
absent. This absence can also easily be rationalized: The 
relaxations of fluctuations in the difference of the mono- 
layer density fields require longer time scales than the 
relaxations of pure shape fluctuations. The slower mode 
no longer describes relaxation of bending fluctuations, al- 
though it can not be approximated simply by ei- 
ther. Clearly, case b) provides a new scenario that is dif- 
ferent from the "conventional" scenario depicted in case 
a). Indeed, a recent analysis of data obtained from flicker- 
noise measurements of giant quasi-spherical vesicles made 
of stearoyl-oleoyl-phosphatidylcholine Q, which seemed 
to have resolved two modes with corresponding time scales 
quantitatively similar to those quoted in Table 2, suggests 
reasons for searching for this scenario. 

We have in the above discussed our interpretation of 
the monolayer compressibility modulus, fc e ff, and proposed 
that the values of this parameter should be lower than 
the standard values quoted in the literature. The numer- 
ical analysis illustrates clearly possible consequences that 
a reduction in the value of k c g can have. As a case of com- 
parison with case b) (Fig. 2) Figure 3 displays the vari- 
ous dispersion relations obtained when the value of k e g is 
increased to 100 erg/ cm 2 and all the rest of the param- 
eters have the same values as in case b). Quantitatively, 
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11525 


4940 


2648 


1583 


1018 


691 


490 


359 


271 


t+ [ms] 


132 


66 


39 


26 


19 


14 


11 


9 


7 


P- 


0.994 


0.987 


0.981 


0.976 


0.970 


0.965 


0.960 


0.954 


0.949 



Table 1. Numerical values of the time scales, r=p = 1/S1^, and the amplitude /3_, characterizing the two dissipative modes 
of vesicle dynamics for «: a = 10~ 12 erg = 24fc B T, rj = 0.01 erg s/ cm 3 , A = 2 x 10~ 5 erg/ cm, E — 10, R — 20 nm, feefi = 
100 erg/ cm 2 , and b = 10 7 ergs/cm 4 . 
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10 


T- [ms] 


31977 


21796 


15485 


11378 


8625 


6720 


5362 


4366 


3617 


r+[ms] 


3185 


995 


448 


241 


146 


95 


65 


47 


35 


P- 


0.287 


0.187 


0.143 


0.117 


0.099 


0.087 


0.077 


0.069 


0.062 



Table 2. Numerical values of the time scales, r=p = 1/S1^, and the amplitude /3_, characterizing the two dissipative modes 
of vesicle dynamics for « a = 10~ 12 erg = 24fcsT, 77 = 0.01 erg s/ cm 3 , A = 2 x 10 -5 erg/ cm, E = 10, Ro = 20 /im, 
keg = 30 erg/ cm 2 , and b = 2 x 10 8 ergs/ cm 4 . 




the increase in the value of fc e ff reduces the characteris- 
tic time scales of both of the modes, affecting the slower 
mode more significantly, though. There is also a qualita- 
tive change in the nature of the two modes for the first few 
low £ values. For the higher value of k c s : the slower mode 
reflects predominantly the relaxation of bending deforma- 
tion governed by the renormalized bending rigidity K a , and 
the faster mode reflects largely the relaxation of the mono- 
layer density-difference field, in reversal to case b). This 
point serves to underline the need to resolve quantitatively 
the issue of whether the renormalization by fluctuations 
of compressibility moduli of lipid-bilayer membranes in 
systems of giant vesicles is significant enough to be ex- 
perimentally relevant. To be sure, reductions in the value 
of fc e ff do not always cause the kind of quantitative as 
well qualitative changes in the relaxation dynamics just 
discussed. In cases where the intermonolayer friction co- 
efficient, b, has lower values, e.g., b = 10 7 ergs/ cm 4 , the 
changes in the relaxation dynamics of the two modes for 



low I values are not significant for the same increase in 
the value of k e g. 

Another point to discuss in connection with experi- 
mental studies of vesicle fluctuations by the technique of 
flicker-noise measurement and analysis concerns the ef- 
fects that variations in the parameter S have on the re- 
laxation dynamics of the fluctuations, as S is a vesicle- 
specific parameter that typically varies in the range from 
to 25 pj] . We have checked systematically in the numer- 
ical analysis the effects, which are summarized by a com- 
parison of Figure 4 with Figure 2. It is easy to see that 
a reduction in S from 10 to leads to a significant in- 
crease in the characteristic time scale of the slower mode, 
but very small increase in the time scale characterizing 
the faster mode. Moreover, the relative positions of the 
different dispersion relations do not change qualitatively. 

For the sake of completeness of the discussion on the 
dissipative dynamics, we mention an analytical relation- 
ship between fi-{t) and which may be exploited 
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Fig. 4. The same dispersion relations as those illustrated in 
Fig. 2, for the same parameter values but one: £ — 0. 



in analysis of experimental data. This relationship and its 
potential have already been pointed out independently in 
Ref. [pi . It follows straightforwardly from Eq. (|9|) that 



£(£ + !)(£ -!){£ + 2) n a 



£(£■ 



21% 



V R 3 r e 
i) fhs 
b 



[£{£ + 1) + S] 



(100) 



In the above expression, all quantities may be considered 
known except the ratio k e g/b. Thus, for vesicle systems 
where the two relaxation processes can both be resolved 
experimentally, the product C2-(£) J? + (£), and in turn, 
fc e ff/6, can be quantitatively determined. This quantitative 
information may be used to obtain a quantitative measure 
of the intermonolayer friction coefficient b, provided that 
fc e ff can be obtained by use of other i ndep endent methods, 
e.g. the micromechanical technique [£5j. This method of 
determining b can serve both as an alternative to, and a 
check of, the method used previously for estimating the 
value of b [^4| . 

Finally, we now return to the issue of the validity of the 
no-inertia approximation, which has already been men- 
tioned in Section 4. This discussion is relevant, especi ally 
in view of the recent publication of Pott and Meleard |l8| , 
in which the experimental data seems to suggest that the 
dynamics of quasi-spherical vesicles is not purely dissi- 
pative, but that the relaxation is modulated by an oscil- 
latory component of a small, but observable amplitude. 
Recapitulating Eq. (M^), and using pt, = 1 g/cm 3 and 
rj = 0.01 erg • sec/cn^for water, we have 

7^10- 6 — . 
to 

Setting a — 60 which corresponds to the lowest spher- 
ical harmonic mode, and using the shortest time scales 



given in Table 1 and Table 2, for the two different cases a) 
and b), respectively, we get two corresponding quantita- 
tive estimates, 1Z ~ 0.04 and TZ ~ 0.001. The issue is then 
really whether these numbers can be effectively treated 
as zero. A theoretical analysis that may address this is- 
sue clearly must take into account the inertial term in the 
bulk hydrodynamics and is outside the scope of this pa- 
per. But, if the approximation turns out to be valid, then 
the theory we have presented here does not provide an 
explanation for the observed oscillatory behaviour. One 
may indeed question whether the observed oscillatory be- 
haviour is genuine. If, however, the oscillatory behaviour 
is genuine, then the no-inertia approximation should be 
re-examined by studying the effect of the inertial term on 
the vesicle dynamics . Moreover, one needs to bear in mind 
that, if it is indeed the inertial term that is responsible for 
the oscillatory behaviour, the time scales that character- 
ize the part of the exponential relaxation of the dynam- 
ics should in principle be qualitatively different from the 
times scales obtained from a theory of purely dissipative 
dynamics such as the one presented in this paper. Looked 
at from this point of view, the interpretation given by the 
authors of Ref. |l8| of their experimental data, which was 
based on their earlier theory of purely dissipative dynam- 
ics fl7|l , does not appear consistent with the same data 
which shows the presence of an oscillation. 



For the dynamics of the density-sum field, whose purely 
dissipative dynamics would be characterized by a time 
scale of the order of 10~ 7 seconds, it is certain that the 
no-inertia approximation does not hold. Consequently, the 
dynamics of this field will contain an oscillatory compo- 
nent. We may argue, however, that this dynamics will 
have no observable effect on the dynamics of the other two 
modes. One argument is that the coupling of the density- 
sum field to the other two fields implied by the free energy 
is very weak for giant vesicles (and is in fact non-existent 
for planar membranes, at least within the framework of 
linearized theories). Another argument is that the time 
scales associated with the relaxation-oscillation dynamics 
of the density-sum field would still be very short compared 
to the time scales characterizing the dynamics of the other 
two fields. 



We would like to end this paper by stating our two 
main hopes that come with the presentation of the work. 
First, the current theoretical work will encourage more 
systematic and careful experimental work that will lead 
to a clear and quantitative understanding of the dynam- 
ics of vesicle shape fluctuations; second, the fact that the 
theory is built upon certain fundamental principles and 
considerations, together with the systematic nature of its 
formulation, will be exploited in terms of extensions of the 
theory that will describe vesicle systems formed of fluid bi- 
layers containing more than one molecular species (e.g. a 
second lipid or a protein). 
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A Appendix 

In this Appendix, we will sketch briefly the calculations 
of static correlation functions based on the free-energy 
model, Eq. (|l2|) . Eq.(p"2|) can be expressed in a form more 
convenient for the calculations presented here: 



JdA {fc cff (4> A ) 2 + k cS (4 

+2\H4> A + a } , 



E f + 2k cS H 2 



(101) 



where the use of the "hat" over the variables indicates 
that they should be considered as stochastic variables. 

The calculations of the static correlation functions are 
essentially based on a Gaussian theory of the fluctuations. 
An inevitable issue in the calculations concerns the choice 
of an appropriate set of independent degrees of freedom, 
whose fluctuations must be integrated over. Although ug rn , 
and i>f ml as defined in Eqs.(fil3|)-(f42|), respectively, 
are independent of each other and appear directly in the 
free-energy expression, they are not appropriate in gen- 
eral. To make the point clear, consider a local area el- 
ement of membrane surface, Aa, associated with which 
are and lipid molecules in the outer and inner 
monolayer, respectively. Changes in the local monolayer 
density fields, <f> A and 4> s can then be brought about by 
two types of independent physical processes: lateral flows 
of lipid molecules at fixed surface geometry, which lead 
to changes in and at fixed A a', and, a change in 
the surface geometry at fixed N\ and , which changes 
A a- Thus, it is more meaningful to choose, together with 
the u — field, the fluctuating fields that reflect the former 
process, as the appropriate set of independent degrees of 
freedom. Such fields, denoted henceforth by n ± (0, </?), are 
related to the apparent density fields as follows: 



p±{9^) = U + pan^e,?)] 



1 



p n ± (9,^)}[l+g 2 {u)} 



j q -r hoi" vy> wjl- 1 t i (102) 

where g(u) = a(u)/R^ and ^(w) = — 2u + 3u 2 — (Vlw) 2 /2 
is the expansion of 1/ \J g{u) — 1 up to second order in u. 
It follows easily that 

4 A (9,<p) = cf> A + n A (0,<p) + <t> A g 2 (u) +n A (9, V )g 2 (u) 
j> S {6,y) = h E (9,<p)+g 2 (u)+n E (9,f)g 2 (u) , (103) 

where n A (9, <p) = (h + (9, if) - n~(9, <p))/2 and h s (9, ip) = 
{n+(9, V )+h-(9, V ))/2. 



Expanding the free energy given in Eq.(lOl) to second 
order in the fluctuations degrees of freedom, eliminating 
WoO) "oo an d ^oo usm g the fixed- volume and the fixed- 
molecular number constraints 



1 1 



\Utr> 



1>1 m= — l 



l A = Const. , 
inn = Const. , 



(104) 

(105) 
(106) 
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leads to the following expression: 
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F = Const. 



{ \ {(* - 1) i l + 2 ) W + l)«-ff + CT oi?o 

- Ks (^) 2 + 4Ai? ^] + 8k s [1 + (<^) 2 ]} \u lm \ 2 
+k s (nfj 2 + Ks (nfj 2 - [{I - 1)(£ + 2)XRo 
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where n s = k c gR,Q. 

Based on the above Gaussian theory for the fluctua- 
tions, and noting that, to the first order in fluctuations 

7 A - A r>lA ~ 

we finally obtain the following two static correlation func 
tions, ^ 



knT 



\Uir, 



(109) 



(»< m %m) = ^ (Nro| ) • (HO) 



which has been quoted in Eq.(|97|) and Eq.(|9S|), and 



4 Other correlation functions can be obtained just as 
straightforwardly, but they are not listed here. 
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